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Abstract We compare results from the traditional partial wave treatment 
of deuteron electro-disintegration with a new approach that uses three di- 
mensional formalism. The new framework for the two-nucleon (2N) system 
using a complete set of isospin - spin states made it possible to construct sim- 
ple implementations that employ a very general operator form of the current 
operator and 2N states. 



O ■ 1 Introduction 

The theoretical description of electromagnetic processes is based, like most 
nuclear physics formalisms, on a partial wave decomposition of relevant oper- 
ators. This restricts theoretical methods to systems where a relatively small 
number of partial waves is important. Recently, several formalisms for the 
three dimensional description of few-body systems and processes have been 
developed. In this paper we fill a gap in the development of the three dimen- 
^{■^ ' sional framework and present an approach that allows for a simple implemen- 
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tation of three dimensional electromagnetic currents. Our final expressions 
can be translated to a numerical implementation via direct substitutions of 
16 dimensional square matrices representing operators in the 2N isospin - 
spin space. Finding the matrix representation of relevant operators is greatly 
simplified by using symbolic programming in the Mathematica ® [I] soft- 
ware package. Our approach allows us to use a very general operator form of 
current operators and can therefore be used for a wide class of processes. 

In this paper we apply it to the case of electron induced deuteron disin- 
tegration and compare the results with traditional partial wave calculations. 
It is worth noting that the methods presented in the following sections can 
be applied to the description of other processes, where electroweak probes 
interact with the 2N system. Electromagnetic currents can be replaced by 
any operators acting on the same degrees of freedom; this makes our imple- 
mentations also useful for calculations involving for example muon capture 
or neutrino induced deuteron disintegration, performed recently still with the 
use partial wave expansion [2H3]- 

2 Formalism and notation 

We adopt a notation in which capital letters describe the total momentum of 
a two particle system (P = pi + P2), lower-case letters describe the relative 
momentum (p = I (pi — P2))- Subscripts denote individual particles and 
superscripts assign a momentum to a particular quantum eigenstate. The 
two particle momentum eigenstates are normalized such that: 

<p'P'|pP)=5 3 (p'-p)5 3 (P'-P), (1) 

J d 3 pd 3 P I pP)(pP |= II (2) 

and the transition from the total and relative to the individual momenta can 
be achieved using: 

Pl = P + ip 

P2 = ~ P, (3) 

where in © and in the following we assume that the difference between the 
proton and neutron mass is negligible. 

We examine the case of deuteron disintegration (e + 2 H — > e + p + n) 
where the 2N system is treated in the non-relativistic approximation. In the 
initial state the deuteron is at rest (P = 0) and the electron has a momentum 
magnitude of q e . We assume that the rest mass of the electron is negligible in 
comparison to its kinetic energy therefore the initial electron energy E e ss q e . 
The final electron momentum magnitude is q' e , the final energy E' e w q' e and 
the electron scattering angle is 9 e The magnitude of the three momentum 
transferred to the 2N system is 

Q = y/q2+q?-2q e q' e cosfl e » y/E* + Ef - 2E e E' e cos 6 e . (4) 
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We work in a reference frame, where the momentum transfer is parallel to z. 
In this frame momentum conservation leads to the expression for the total 
momentum of the proton and the neutron in the final state: 



The magnitude of the final relative momentum can be calculated from the 
energy conservation: 



where E d is the (negative) deuteron binding energy and the direction of p f 
can be arbitrary. The crucial nuclear matrix element 1VP between the initial 
deuteron state (where the total momentum P = and the two particle total 
angular momentum has a z projection wid) and the final 2N scattering state 
can be expressed in terms of the full 2N current operator (j 2Ar ) and the t 
operator: 



M" (p f ,P f ) 

= a (p f P f ,mifi,m2^ | {% + t(E)G (E))j% N \ q> d m d P = 0) 
= Q (p f P f ,mi^i,m2^ | {l + t{E)G Q (E)) 
(f(l)+f{2)+f{l,2)) \(f> d m d P = 0) 
= 2 a (p f P i ,m 1 v 1 ,m 2 is 2 |^(2) | (f) d m d P = 0) 
+ a (p i P\m 1 v u m 2 v 2 \.f(l,2) \(f) d m d P = 0) 
+ a (p f P f ,TOi^i,m 2 ^ 2 | tG j 2N I 4>dm d P = 0). (7) 



In this equation the final state is anti-symmetrized: 

a (pP,miy 1 ,m 2 !/ 2 |= - ((pP, mivx, m 2 v 2 | -(-pP, m 2 i'2, "H^i I)- (8) 

In the first term on the right hand side of © particle 1 has a spin (isospin) 
z projection mi {vi), particle 2 has a spin (isospin) z projection m 2 (v 2 ) and 
in the second term on the right hand side of © the particles are exchanged. 
The j(l) 0(2)) operator is a single nucleon current acting on the degrees 
of freedom of the first (second) nucleon. The 2) operator accounts for 
processes where two nucleons are involved, t is the 2N transition operator 
and Go is the free 2N propagator. The energy argument of the transition 
operator and the propagator is E = (p l ) 2 /m. Finally, the \i index denotes 
the component of the current operator. In particular [i = stands for charge 
density operators while /i = 1,2,3 stand for spatial components. In the fol- 
lowing sections the way we calculate the individual parts of M M in ([7]) will be 
discussed separately for a specific choice of the coordinate system and value 
of /i; we will drop /i for brevity. 



P f = Pi +p 2 = (0,0,Q). 



(5) 



|p f I - I^Pi - P2)l = -^Am{E d + E e ~E' e )~Q\ 



(6) 
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3 Deuteron bound state 

The structure of the deuteron wave function can be written in the operator 
form, following (HEJE)]: 

| (f) d m d P = 0) = 

2 

d 3 p | pP = 0) £ <MIP|) [I™ ® &i(p) spin ] [| 00) ® X (m d )} 
i=i 

2 

d 3 p | pP = 0) M\P\) [BiiP)} 00) ® X (m d )] . (9) 



i=i 



where 



bi = 1 

<r(l) • p<r(2) • p - ip • pi 

In equation Q | x( m d)) i s a state in which the spins of the two spin ^ 
particles are coupled to a total spin 1 with a z projection m^. The isospins 
of the two nucleons are coupled to the total isospin state | 0). <fii are 
scalar functions of the relative momentum and <x(l), <x(2) are doubled spin 
operators in the spin space of one nucleon and identity operators in the spin 
space of the other nucleon, respectively: 

<r(l) = (ff i 8l,fl*0l,ff , 8l) 1 (10) 

er(2) = (I® o^lBo*,!®**). (11) 

Vector components in (fTCTI) and (ITT]) are operators written in terms of the 
tensor product <£>. They act in the 4 dimensional spin space of the two nucleon 
system and can be represented by 4 x 4 matrices - tensor products of identity 
operators and Pauli matrices. 

Note that in equation © states and operators in the 2N isospin - spin 
space are placed inside the [. . . ] brackets; in Sec. [7] we will introduce a simple 
way to implement these expressions as 16 dimensional vectors and 16 x 16 
matrices for use in our numerical treatment. 

Scalar functions <fii in expansion Q can be calculated using three di- 
mensional formalism, see for example [B]. Nowadays deuteron bound state 
calculations can use any 2N potentials given in the operator form and do not 
require substantial computational resources. 



4 Single nucleon currents in three dimensions. 

Single nucleon (IN) currents act on the degrees of freedom of one particle. 
Their matrix elements in the momentum space depend only on the initial 
and final momenta and are operators in the isospin - spin space. For example 
the matrix element for the second nucleon, j(2), reads: 

[(p'iP 2 |i(2) |pipa)] =a(p'i-Pi)b'(2,P2-P2,P2+P2)] (12) 
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where in view of the standard nonrelativistic current, the dependence on the 
difference and sum of the initial and final momenta is used. Implementing 
the transition from the individual particle momenta to the relative momenta 
leads to 



[<p'P'| j(2)|pP)]=5(V-ip + p'-p 



p) 



(13) 



Again, the expressions inside the square brackets [. . . ] can be easily repre- 
sented using the notion of the Kronecker product, see Sec. [7J The action of 
j{2) on the deuteron state at rest can be worked out: 



[(p'P' I j(2) \<j> d m d V = Q)] = 

= E<MIP' + ? P 'IMJ(2,P',-2P')] 



i=i 
— IY-|1N 



W+2 P ') 



[| oo)®xKz)] 



[0 1N (2,p',P')] [| 0>®x(m d )] 



(14) 



where we used P = and the normalization of momentum eigenstates (JT])-([2]) 
so 1N is the resulting single particle operator. Equation ((T4)) gives the full 
isospin - spin state for the final p', P' momenta. 



5 2N currents in three dimensions. 

For a wide class of 2N current operators, their matrix elements in the mo- 
mentum space (operators in isospin - spin space) are given in the form: 

[(p' lP ' 2 | j(l, 2) | p lP2 >] = [7(1, 2, p[ - Pl , p' 2 - pa)] , (15) 

see for example [7JE1IHI ■ The right hand side of (TT5|) is as a linear combination 
of scalar functions (// , ff) and products of spin space operators (Ojs , 
Oj) and isospin space operators (Tj): 

[j°(l,2)] =J2Y,fi 3 (Pi -P1.P2- Pa) [TiO iS ]. (16) 

i=l j=l 
5 24 

[j(l,2)] =^^/f(pi - P1)P2 - pa) [T % 3 ] , (17) 

i=l 3=1 

where the subscript S distinguishes density and current components. In Ref. 
[ID] a general operator basis for the local 2N current operator was intro- 
duced. In Sec. [JJ some examples of operators from (fT6|) . ([17]) will be used to 
demonstrate our matrix representation of expressions inside [...]. 
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Again, using ([3]), current matrix elements become: 
Kp'P'l j(l,2)| P P)] 



j(l,2,ip'-ip + p'- P) Ip'-Ip-p' + p) 



( 18 ) 

In the following we restrict ourselves to this class of momentum dependences. 
Our approach can, however, be generalized to include any type of momentum 
dependence. The action of 2) on the deuteron state can be worked out 
using ((I]),©, © and P = 0. In the laboratory frame it yields: 



dV5>(i P "l) 



7 = 1 



Kp'P'|j(l,2)|0 d m d P = O)] 
j(l,2,ip' + p'-p",~P'-p' + p'' 

[B,(p")][|0Q)®x(m < i)] 



= [0 2N (1,2, P ',P')] [|00)» X W], (19) 

where (9 2N is the resulting two-particle operator. Equation (fT9")l gives the full 
isospin - spin state for the final p', P' momenta. 



6 t operator in three dimensions 

The t operator satisfies the Lippmann - Schwinger equation: 

t(E) = V + t(E)G (E)V (20) 

or equivalently 

t{E) = V + VG {E)t(E) (21) 

where Gq(E) is the free propagator depending on the energy E and V is a 
2N potential. It follows that, as shown in [11] also t can be written as a linear 
combination of scalar functions (t 7 ,i) and operators (W 7i j) in the isospin - 
spin space: 

[<p' | t(E) | p>] = J2 E hAlv'l IPl, P' ■ P, E) [W 7li (p', p)] . (22) 

7 i= 1 



Here 



W 7>i ( P ',p)] = C~®wf n (p',p) 



are again operators in the isospin - spin space (matrix elements between 
momentum states), with w, (i — 1, 2, 6) acting in the 4 dimensional spin 
space of the 2N system. The decomposition (|22l) is not unique; our choice 
of the six Wi operators is consistent with [6] . Scalar functions arising in 
the decomposition of t can be calculated in the three dimensional formalism. 
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Calculations can be performed for any type of the NN potential satisfying a 
similar expansion (l2"2"j) . For details see Ref. [5]. 



wi(p',p) = 1 
w 2 (p',p) = <r(l)-<r(2) 
w 3 (p',p) = z(<T(l) + <T(2)).(pxp') 
w 4 (p', p) = o-(l) • (p x p') <r(2) • (p x p') 
w 5 (p', p) = <t(1) • (p' + p) <r(2) • (p' + p). 
w 6 (p', p) = o-(l) ■ (p' - p) <x(2) ■ (p' - p) 



(23) 
(24) 
(25) 
(26) 
(27) 
(28) 



The C 7 =| 7) (7 I isospin operators project onto one of the four 2N isospin 
states: 

liy = Q } lm t = -t...t). (29) 



7) 



The I 7) states are chosen in this way, because t(E) conserves the total 2N 
isospin. 

The rescattering part of the matrix element M in ([7]) can be written as: 

[(p'P' \t(E)G j 2 N |^m d P = 0)] 

1 



J d 3 p[(p'|*(£)|p) 



x[0(p,P')] [\00)®x(™ d )] 

„/9 „9 fl IPI 



Ip'I [f(|p'D] 



In 



P 



P - IP' 

x [| 00)(g>x(m d )} 



(30) 



where O is either 1N from Eq. flU) or 2N from Eq. ([15]). E = p' 2 /m is 
the relative energy of the final 2N state and 



[fdPl 



2tt 



d<t> dx[(p'\t(E)\p)]{0(p,P')], 



(31) 



since 



p = |p| (vl — x 2 cos<p, y/l — x 2 sin0, x). 



(32) 



The integral (|3T)]) with the cut-off value p can be easily calculated numer- 
ically. In the next section we show how to prepare its component [0(p, P')]. 



7 Representation of spin-isospin operators 

Now that the form of expressions in (JTJ) has been established it remains to 
find a way to represent operators and states inside [...]. Once an appropri- 
ate matrix representation is found, numerical calculations can be constructed 
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Table 1 Reference quantum numbers for our KP states. 



1 

2 
3 
4 
5 
6 
7 
8 
9 
10 
11 
12 
12 
14 
15 
16 



2 
1 
2 
1 
2 
1 
2 
1 
2 
1 
2 
1 
2 
1 
2 
_ 1 
2 

_ 1 
2 

_ 1 
2 

_ 1 
2 

_J_ 

2 

_J_ 

2 

_J_ 

2 

_J_ 

2 



^(i) mg°(j) 



» in (i) 



_ 1 

2 

1 

2 

_ 1 

2 

1 

2 
_ 1 
2 

1 
2 
_ J_ 

2 

1 
2 
1 

" 2 ' 



using equations (jT4"]) . (|19l) and ([50)1 by simple substitutions and matrix mul- 
tiplications. 

Our choice for the 16 dimensional basis of the two nucleon isospin - spin 
state space (the deuteron in the initial state, the proton and the neutron in 
the final state) is the set of vector states {| a)} (i = 1, 2, 16): 

| ei) = (| mf spin (*)>® | mr spin M)) ® (| mf in (i))® | mf "(*)>) (33) 

where TO ^P m ( lsos P in ) are ^j^g S pi n (isospin) projections of nucleon j and the 
corresponding quantum numbers are given in Tab. [1] 



Any operator or state in this space can be constructed using the notion of 
the Kronecker product (KP). The Mathematica ® pQ symbolic programming 
software contains definitions for the KP, which makes translating expressions 
for isospin - spin operators a simple task. However, it is important to remem- 
ber that the order of operators in the KP must be preserved. Tab. [T]can serve 
as a reference to keep consistence with this paper. 

The deuteron operators [-Bi(p)] and [^(p)] with P = {Px,P v ,Pz) from 
equation (O have a form simple enough to have their matrix representation 
written out in full. [B\] is simply a 16 x 16 identity matrix. [B2] has a block 
diagonal form: 

/BO \ 

B 

6 
V BJ 
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where B is a 4 x 4 matrix: 

( i (~ p l - p l ~ p l) + p i p * ( px ~ ip y) 

l i ■ \ 1 ( 2 2 2\ 2 

(pa + *Py) § l^-Pa; ~ p y ~ p z) - "P z 

Pz (Px + ipy) {Px — ipy) (px + ipy) 

\ (Px + iPy) Pz (- (Ps + iPy)) 



Pz {Px — iPy) (px - iPy) 2 N 

{Px — ipy) (Px + ipy) Pz (— {Px ~ ipy)) 

i (-Pl - Pl ~ Pl) - Pl P* (- (Px - ipy)) 

Pz (- (P. +ipy)) i (-Pl - P y - P\) + P\ ) 



Operators appearing in (TTB")) and (fTTf have a more complicated isospin 
operator form. A simple example of their matrix representation is given below 
in Fig. [U 
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Fig. 1 


[T 3 Og] (pi - pi 


= (qix,qi y ,qiz) 


,P ; 2 


- P2 


= («2as, 


qiy, qiz)) with T 2 = 


= Tf- 



T2 and O5 = (qi x cri) + (q2 x (T2). t< (<Xi) is the 2N isospin (spin) vector operator 
acting in the space of nucleon i. 



States can be constructed in a similar manner using the built in Mathe- 
matica ® definitions for Clebsch-Gordan coefficients, for example [| 0)® x(l)] 
results in: 

'0, 0, 0, 0,-^,0, 0, 0, - i 0, 0, 0, 0, 0, 0, 0) . (34) 



8 Results 

In the following we will present results obtained using a chiral NNLO po- 
tential [E] with A = 550 MeV/c and A = 600 MeV/c. The operator form 
of such a potential was briefly described in Appendix C of Ref. [5], where 
also an example set of necessary parameters was given for its neutron-proton 
version. The same parameters will also be used in the present paper. 

There are three basic ingredients in our calculations: the deuteron wave 
function, the 2N t-matrix and the 2N current operator. Before we show se- 
lected observables for the 2 H(e, e'p)n reaction, we will describe our numerical 
performance and the way we verify the quality of our calculations. 

As described in Ref. [5] and equation ©, the deuteron in the operator 
form is represented by two functions 4>i(p) an d <fo(p)- The corresponding 
Schrodinger equation for 4>i(p) and <p2 (p) can be rewritten as an eigenvalue 
problem, which is of the same type and dimension as the one solved for the 
deuteron wave function in the standard partial wave representation, where 
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one deals with the s- and d-components, ipo(p) and ip2(p)- The connection 
between the solutions is very simple [5] 

Mp) = V4^Mp), Mp) = ^-^Mp) (35) 

and can be used to check the numerical performance. 

In Fig. [5] we show directly the <pi(p) and 4>2{p) functions required for the 
operator expansion of the deuteron. In Fig. [3] the resulting s- and d-wave 
components in momentum space are compared to the results obtained by 
firstly decomposing the NN potential into partial waves and then solving the 
resulting Schrodinger equation in its standard form. The agreement for the 
two wave function components is perfect for all their significant values. 




0.5 1 1.5 2 0.5 1 1.5 2 

p[fnf'] p[fm- 1 ] 

Fig. 2 The 4>i(p) (left) and 4>2 (p) (right) expansion function in the operator form 
of the deuteron as a function of the magnitude of the relative momentum p for the 
considered chiral NNLO potential. 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 
p[fm" 1 ] 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 
P [fm" 1 l 



Fig. 3 The s-wave (left) and d-wave (right) component of the deuteron wave func- 
tion as a function of the magnitude of the relative momentum p for the considered 
chiral NNLO potential. Crosses represent results obtained using the operator ap- 
proach and solid lines are directly from the standard partial wave decomposition. 



In Ref. ;0j we solved the Lippmann-Schwinger equation (LSE) for the 2N 
t-matrix directly in three dimensions. At that time we focused mainly on 
the on-shell behavior of the expansion coefficients ti(p' ' ,p 1 x; E2n), that is we 

were interested in ti [po,po, x; E2N = m )' which are sufficient to calculate 
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the Wolfenstein parameters and the nucleon-nucleon scattering observablcs. 
Furthermore, we solved the LSE in such a form (Eq. (2.6) from Ref. [5]) 



3 

J d3 P" E V T (P'> P")G (p'Yp (p", p)B kjj , (p', p", p) 



(36) 



that the magnitude of the initial p momentum could be fixed. 

Clearly, for the 2 H(e, e'p)n reaction we need a "left" version of Eq. (|3"6"1) . 
which allows us to find the half-shell t-matrix for a fixed final relative momen- 
tum, p f , given now by the reaction kinematics. The starting point for this new 
version is equation (|2"0|) . Repeating the algebra outlined in Ref. [B], we pre- 
pared a numerical realization of this "left" version of the LSE, leading to the 

scalar expansion coefficients ^ ^p f ,p; E 2 n = ) = t% (p 1 \p \x; E2N = 

where x = p -p. Our numerical scheme was again based on matrix inversion 
and used the standard LU decomposition of Numerical Recipes [T3] • In order 
to achieve a unique and smooth solution also for p = p { , it was sufficient to 
calculate the average 

( (P f ) 2 

U p ,p,x\E 2 N - 



in 



If ( (P f ) 2 \ f (P f ) 2 

= ^ ( M P l ~ V ' P) x ! E in = J + ( P l + S p f , p, x; E 2 n = 

(37) 

with 8 pt « 0.01 fin" 1 . 

Actually, this effort turned out to be unnecessary and provided merely 
an additional check of numerics, since 

U (p l ,p, x; E 2N = = U (p,p l , x- E 2N = (38) 

for the most general rotational, parity and time reversal invariant form of 

the NN force. That means that the left coefficients, U (pf ,p, x; E 2 n = ^ , 

can be obtained directly from the "right" version of LSE. 

In order to further check our t-matrix coefficients, we used them to cal- 
culate the partial wave representation of the t-matrix: 

J\2 - 



tf(l's)j I t (e 2N = b±) \p(ls)j), 



where I (l 1 ) is the initial (final) angular momentum of the 2N system, s is the 
2N (conserved) spin and j is the 2N (conserved) total angular momentum. 
These matrix elements can be compared with direct solutions of LSE obtained 
in the standard partial wave representation. We performed the projection 
of the three-dimensional t-matrix on partial waves, employing the simple 
method proposed in Ref. [14] for NN forces. In Figs. |4][7]we show examples 
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for the uncoupled and coupled channels, with the 2N isospin t — and 
t = 1. We chose p ~ 1.9 fm -1 , which corresponds to a relatively high NN 
center of mass energy, E2N— 150 MeV. For such an energy many partial 
waves contribute to the NN scattering observables and the question arises 
if the partial contributions are consistent with the full three-dimensional 
calculations. From Figs. HE (and many other cases which are not shown 
here) we infer that this is really the case. The agreement between results 
based on the two quite different approaches is very good. 





Fig. 4 The real (left) and imaginary (right) parts of the half-shell 1 So t-matrix as 
a function of the initial momentum p for p l ~ 1.9 fm _1 . Points represent predictions 
obtained by a projection from the three-dimensional results. Solid lines represent 
direct solutions of LSE in the standard partial wave decomposition. 





Fig. 5 The same as in Fig. U for the half-shell 3 H$ t-matrix. 



The final ingredient in our framework is the 2N current operator. It con- 
sists of the single- nucleon and two-nucleon operators. For the purpose of this 
paper we assume that its single-nucleon part comprises the standard nonrcl- 
ativistic charge density as well as the convection and spin current operators. 
In the 2N part we take for simplicity only the leading one-pion-exchange cur- 
rent operator in the chiral effective field theory representation. In our three- 
dimensional treatment of the 2 H(e, e'p)n reaction, we calculate the spin and 
isospin matrix elements of the current operator directly, using simple ma- 
trix representations of the spin and isospin operators and the concept of 
the Kronecker product to deal with the two-nucleon spin and isospin spaces. 
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In the traditional calculations, a partial wave decomposition of the current 
operator is required. It is a rather easy task for the single-nucleon part of 
the current. For the one-pion-exchange current operator it is known analyt- 
ically (see for instance Ref. [H]). It can also be obtained using the method 
proposed in Refs. |16U17j . where even more complicated two-pion-exchange 
current operators were considered. 

To give examples of our results on deuteron electro-disintegration we 
chose several electron kinematics given in Table [5J They allowed us to study 
the reaction for three different internal nucleon-nucleon energies (correspond- 
ing to the three values of p i ) and for five values of the three- momentum 
transfer Q. The first parameter is the input for the t-matrix calculations and 
the second one specifies the properties of the current matrix elements. 

Table 2 The six electron kinematics considered in the paper for the exclusive 
2 H(e,e'p)n process. The initial electron energy (E e ), the electron scattering angle 
(0e), the final electron energy (E' e ), the final relative nucleon-nucleon momentum 
(|p f |), the energy transfer (lj), and the magnitude of the three-momentum transfer 
(Q) are given. 





E e 


8e 


K 


p* of Eq. © 


w = E e - E' e 


Q of Eq. @ 




MeV 


deg 


MeV 


MeV/c 


MeV 


MeV/c 


AT 


500 


6.9 


490.3 


78.1 


9.7 


60 


K2 


500 


17.4 


485.3 


78.1 


14.7 


150 


K3 


500 


6.1 


467.0 


158.7 


30.0 


60 


K4 


500 


36.4 


447.0 


158.7 


53.0 


300 


K5 


500 


16.3 


337.1 


375.3 


162.9 


200 


K6 


500 


73.7 


281.2 


375.3 


218.8 


500 



In addition to the information given in Table [2j we need to label the 
exclusive kinematics. For the fixed "electron arm", we deal in fact with a 
two-body kinematics in the final proton-neutron system. We restrict ourselves 
here to the case where the proton is ejected in the electron plane, where 9 P 
would be the angle between the three-momentum transfer Q and the final 
proton momentum p p . Since we have to distinguish between the (j) p =0 deg 
and (/> p =180 deg cases, we ascribe the negative sign to 9 P for </> p =0 deg. This 
is shown in Fig. [5] Note that the six electron kinematics provide a unique 
solution for any 9 p value and that 9 p changes from do 180 degrees. 

We are now ready to show our results for several selected observables. We 
chose first of all the unpolarized cross section, d 5 a/(dE e > df2 e > df2 p ). We take 
also into account one example of the spin-dependent helicity asymmetry, 

= ojh = +l,3 d )-a(h = -1, J d ) 
11 - a(h = +l,J d )+a(h = -l,J d ) 

where h is the initial electron helicity and the projection of the initial deuteron 
total angular momentum (J^) on Q, Jd z , is equal 1. In addition we show our 
predictions for the deuteron tensor analyzing powers Tk q . Note that they are 
calculated in the system, where Q || z. 

Our primary goal was to compare results based on the partial wave de- 
composition for the t-matrix and the nuclear current operator with new pre- 
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dictions resulting from the three-dimensional scheme. We observed a perfect 
agreement for all the electron kinematics and for all the considered observ- 
ables, if a sufficient number of partial waves in the first type of calculations 
is included. The six kinematics can be divided into two groups: (Kl , K3 , 
K5) and (K2 , KA , K6). In each group a similar type of convergence of the 
observables with respect to the number of partial waves is observed. That is 
why in Figs. IMT^l we show predictions for two representative (K3 and K6) 
kinematics only. In the first case we see a rapid convergence and partial wave 
based results with j < 4 are already very close to the full three-dimensional 
prediction. In the second case all partial waves with j < 9 are necessary to 
achieve convergence. 

It is interesting to see that slow convergence for the K2, K4 and K6 kine- 
matics does not result from the higher p f values (that is from the t-matrix) 
but is related to the Q values and thus to the partial wave decomposition 
of the current operator. It is well known (see for example Ref. [TB]) that 
especially the partial wave decomposition of the single nucleon current re- 
quires many partial waves. However, even if the initial bound state is given 
in the partial wave representation, the single nucleon current can be applied 
directly in the case of the plane wave amplitudes. This holds not only for 
the two- but also for the three- nucleon system [T5]. In order to demonstrate 
this behaviour, we showed in Fig. [15] observables for the K6 kinematics. In 
this case the single-nucleon current contribution to the plane wave ampli- 
tude is calculated without partial wave decomposition. We clearly see that 
the convergence is significantly improved, even if the two-nucleon current 
contribution to the plane wave part of the nuclear matrix element and the 
whole rescattering part of the nuclear matrix element is calculated with the 
partial wave decomposition. 
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Fig. 6 The same as in Fig. [3] for the half-shell 3 P2 — 3 Fb t-matrix. Rows show 
different I and l' cases (from top to bottom): (7 = 1, l' = 1), (I = 3, I' = 1), (I = 1, 
I' = 3) and (7 = 3, I' = 3). 



16 





5e-07 


-5e-07 

-1e-06 
-1.5e-06 

-2e-06 
-2.5e-06 

-3e-06 
-3.5e-06 

-4e-06 
-4.5e-06 





1e 


06 r 






« 










-1e 


06 - 




-2e 


06 - 


CL 








-3e 


06 - 




-4e 


06 - 




-5e 


06 - 










-6e 


06 - 


E 


-7e 


06 - 




-8e 


06 - 



1 2 
Plfrn 1 ] 



Fig. 7 The same as in Fig. U for the half-shell 3 He — 3 Je t-matrix. Rows show 
different I and I' cases (from top to bottom): (I = 5, l' = 5), (7 = 7, Z' = 5), (7 = 5, 
I' = 7) and (Z = 7,1' = 7). 
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Fig. 8 The kinematics for the exclusive 2 H(e,e'p)n process, k (k') is the initial 
(final) electron momentum. We neglect the electron mass, so | k |= E and | k' | = 
E'. 
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Fig. 9 The unpolarized cross section d 5 a/(dE e t dQ e i dfi p ) as a function of the 
proton scattering angle 9 P for the K3 electron kinematics from Table [2] In the 
left panel plane wave results (dashed line) are compared with results of the full 
calculations (solid line) obtained within the same three-dimensional scheme. In the 
right panel convergence of the full results calculated with a different number of 
nucleon-nucleon partial waves towards the full three-dimensional prediction (solid 
line) is shown. Partial wave based results with j < 1 (dash-dotted line), j < 2 
(dotted line) and j < 4 (dashed line) are displayed. 
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Fig. 11 The same as in Fig. [5] for the deuteron analyzing powers T^ q . 
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Fig. 12 The same as in Fig. for the K6 electron kinematics from Table [5] In 
the right panel partial wave based results with j < 4 (dash-dotted line), j < 7 
(dotted line) and j < 9 (dashed line) are compared with the full three-dimensional 
prediction (solid line). 




Fig. 13 The same as in Fig. [T^Jfor the spin-dependent helicity asymmetry A\\. 




Fig. 14 The same as in Fig. [T^Jfor the deuteron analyzing powers Tk q . 
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Fig. 15 Various observables for the K6 electron kinematics from Table [2] The 
single- nucleon current contribution to the plane wave part of the nuclear matrix el- 
ement is now calculated without partial wave decomposition. The remaining parts 
of the nuclear matrix element are still calculated using partial waves. Results ob- 
tained with j < 1 (dash-dotted line), j < 2 (dotted line) and j < 3 (dashed line) 
are compared with the full three-dimensional prediction (solid line). 



9 Conclusions and outlook 

The presented method to treat several electroweak processes involving IN and 
2N current operators in three dimensions can successfully replace standard 
partial wave treatment. We showed, for the case of electron induced deuteron 
disintegration, that results obtained using the new approach agree very well 
with those obtained using PWD. For all observables considered in this paper, 
the traditional results converge to the three dimensional predictions when the 
number of partial waves is sufficiently high (Figs. I MT1)) . 

Our formalism employs the two nucleon bound state, the 2N t matrix and 
the current operators in the joined isospin - spin space of the 2N system using 
the three dimensional formalism in 2N momentum space. Each element of 
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this framework has been separately tested and compared with the standard 
PWD approach. Our method seems to be more flexible and can deal with the 
rich structures of the 2N force and the current operator, especially derived 
within the higher orders of the chiral effective field theory [T!?lrTWT(T] . We plan 
to use our framework for other processes such as muon capture or neutrino 
induced deuteron disintegration. In [2 muon capture on 2 H and 3 He was 
considered using the PWD approach. It would be interesting to compare 
those results with three dimensional calculations. A similar convergence to 
three dimensional results as in Figs. l9HT4lis expected. 
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